Shape of Clusters as a Probe of Screening Mechanisms in Modified Gravity 
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Scalar fields are crucial components in high energy physics and extensions of General Relativity. 
The fact they are not observed in the solar system may be due to a mechanism which screens their 
presence in high dense regions. We show how observations of the ellipticity of galaxy clusters can 
discriminate between models with and without scalar fields and even between different screening 
mechanisms. Using nowadays X-ray observations we put novel constraints on the different models. 



General Relativity (GR) is a successfully tested theory 
in solar system scales and bellow. The assumption that 
this theory is valid also at cosmological scales gave rise 
to the ACDM model, which has its foundations in two 
unknown components: dark matter and dark energy. The 
nature of these two components could therefore be an 
indication of the breaking down of Einstein's gravity on 
large scales. This has motivated the proposal of several 
theories which modify GR at astrophysical scales 

An imperative requirement to all Modified Gravity 
proposals is that they all must recover GR in the so- 
lar system. This is done via a screening mechanism. 
Presently, there are three main screening mechanisms: 
Vainshtein Q, Symmetron Q and Chameleon In 
this letter we focus on the two later ones, since they are 
described by scalar degrees of freedom, have the com- 
mon feature of emerging at the onset of nonlinear struc- 
ture formation. Moreover, they are used in a wide range 
of models which extend Einstein gravity: f(R), scalar- 
tensor, vector-tensor theories, etc. 

The key feature of screening mechanism is to switch off 
the extra degrees of freedom inside matter overdensites 
(small scales) , and to switch them on in the cosmological 
background (very large scales). When the scalar fields 
are on a fifth force emerges between the matter parti- 
cles. When it is off (the field is screened) the fifth force 
disappears and GR is recovered. This is a highly nonlin- 
ear process, since the scalar fields are strongly coupled to 
matter and have highly nonlinear bare potentials. 

The aim of this Letter is to investigate signatures that 
the chameleon and the symmetron fields imprint in the 
formed nonlinear structures, which can be measured by 
nowadays' experiments, and so be used as probes of ex- 
tensions of GR and the required screening mechanism to 
reproduce it within the solar system. 

Dark matter halos are not spherical, and its density 
scales differently along the x, y or z directions. Such 
anisotropic shape of the halos (ellipticity), leads to an 
anisotropic screening mechanism of the scalar fields. The 
result is that the fifth force between dark matter parti- 
cles can be present in one direction, while being almost 
nonexistent in another. This anisotropy may lead to an 
increase in the ellipticity of the clusters, which can be 
measured by lensing or X-ray observations, and so be 
used as tracers of the inherent screening mechanism. 



It is well known that the Newtonian potential of tri- 
axial systems acquires a shape that is more spherical 
than the matter-density distribution itself. In the case 
of strongly coupled scalar fields their iso-surfaces are ex- 
pected to follow closely the matter-iso-density contours. 
Due to the screening, the fifth force range and couplings 
change along the matter-iso-density contours, leading to 
modifications in the shape of the system. In this Letter 
we test this conjecture by studying the 3D distribution 
of scalar fields that correspond to triaxial dark matter 
halos, and use X-ray observations to put bounds on the 
models with and without scalar fields and to distinguish 
between different screening mechanisms. The Newtonian 
potential 4>n is given by: 
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where S is the over-density denned as Sp/pb, Pb is the 
mean density of the universe, 4>n is the perturbation in 
the metric, fl m is the mean density of the universe in 
terms of the critical density, Hq is the Hubble constant 
and a is the expansion factor. 

The symmetron model Q is defined by the following 
effective potential: 
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which leads to the following equation of motion in the 
static limit: 
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where p and M are mass scales, p is the matter density 
and A is a length scale. We normalize the field <j> 8 with 
the minimum of the potential <fis,a that corresponds to 
zero density and is given by: </> 2 = By dividing 
the whole equation by (f> s o, defining the dimensionless 
quantity \s 

M 2 p? at z SSB , we get: 



, and taking into account that ,055 g = 
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where rj is the matter density field normalized with the 
mean density of the universe, and X s ,o = -J^- is the range 
for the field that corresponds to zero density. 
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The associated effective potential for the chameleon is: 

V c>af f (</> c ) = Mt +n r c n + pe^/ M - , (5) 

where M c has units of mass, j3 is dimensionless, n is a 
positive constant and M p i is the Planck mass. The lin- 
earized equation of motion for the scalar field is: 
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As in the symmetron model and for numerical conve- 
nience, we normalize <J) C with the minimum of the effec- 
tive potential 4> c ,o- In the chameleon case, the minimum 
diverges when the density goes to zero. Thus, in this case 
we normalize with the minimum that corresponds to the 
mean density of the universe. After including the range 
of the field for this particular density 
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and re-normalizing the field Xc — 4>c/4>c,o, we get 
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where Xc = is chosen to facilitate the numerics. 

To solve the field equations we use a Fourier based 
method for the Newtonian case and an implicit multigrid 
non-linear solver for both scalar field equations. The code 
uses a uniform grid and is an extension of [5] , to which we 
added both scalar field solvers. The boundary conditions 
are periodic. The three solvers were tested successfully 
against analytic solutions for a sphere of uniform density 
located in the center of the box. 

To test the possibility that the presence of a scalar 
field can have an impact in the shape of clusters, we 
calculated the Newtonian potential and both scalar fields 
for a density distribution given by a NFW profile @. The 
virial radius R v of the object was fixed to 1 Mpc/h, which 
corresponds to an object of 10 14 M Q (i.e. a cluster of 
galaxies). Following mass-concentration relations coming 
from simulations [7J, we choose a concentration of 6.3 for 
our halo. The density distribution was defined in every 
node of the grid following the analytic profile. The size 
of the box is 16 Mpc/h and 512 nodes per dimension 
were employed for the calculations, which corresponds to 
a spatial resolution of about 30 kpc/h. The underlying 
cosmology needed to normalize the density profile was 
chosen as ACDM, defined by f2 m = 0.3 and Hq = 70 
km/sec/Mpc. All our analysis are at redshift z — 0. 

As we want to measure how closely the scalar field 
follows the triaxiality of the density distribution, we need 
to assume a density profile. We choose the NFW profile: 
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FIG. 1: Contours of density distribution, Newtonian potential 
and symmetron and chameleon fields for the triaxial halo in 
the plane x — y. The contours where chosen such that all 
of them pass through the point (R v , 0,0). The scalar field 
models are defined by zssb = 1.6 and A Sj o = A c = 1.1. 



but using an ellipsoidal radius 



(10) 



instead of the radius r that corresponds to spherical co- 
ordinates. The axial ratios of the density distribution 
where fixed to be (<? 2 ,s 2 ) = (0.5,0.3). 

Our calculations show that the chameleon model is sen- 
sitive to the behavior of the density far from the center of 
the halo. In order to get stable results, we immersed the 
halo in a background with a constant density of 0.4 the 
mean density of the universe. We find that the Newto- 
nian and symmetron values are independent of the pres- 
ence of this background. 

Given the solutions of the field equations, we now mea- 
sure the shape of their iso-surfaces: the axial ratios of 
the iso-surfaces of a given distribution / (in our case the 
Newtonian potential and both scalar fields) is: 
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As the absolute value of the scalar fields do not decrease 
with radius, we defined / as a renormalized version of the 
fields. These iso-surfaces can be approximated by ellip- 
soids, described by the radial coordinate k (eq. (|10p ) with 



axial ratios q 2 = ^ and s 2 = where M xx , M, 



and M zz are the eigenvalues of My. As our system is 
aligned with respect to the box, we do not need to diago- 
nalize My, since eq. ([Tl"T) already provides the eigenvalues. 
The integral in that equation is computed by making a 
summation over the grid up to twice the virial radius. 
As in the symmetron model the scalar field can be com- 
pletely screened up to values of k of the order of 2R V , we 
extended the domain of the integral for the most extreme 
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FIG. 2: Axial ratios q 2 (left) and s 2 (right) for the Newtonian potential and the two implementations of the scalar field as a 
function of the range of the scalar field. The two curves for the symmetron model correspond to two different values of zssb- 
In the case of the chameleon, we show results for two different values of n. The continuous black lines are the values that 
correspond to the density distribution and the dotted black lines the ones that correspond to the Newtonian potential. 



cases (with very small values of A s .o and zssb)- Thus, in 
the symmetron case, we integrated up to the maximum 
between 2R V and the radius k at which / reaches 0.5. In 
all the cases, the summation was made up to no more 
than 3R V . The shape of the region in which the integral 
is calculated was determined iteratively as in [8J. 

Fig[T] shows an example of the phenomenology we 
are dealing with. The plot shows iso-densities and iso- 
surfaces for the Newtonian potential and both scalar 
fields in the plane that corresponds to the major and 
intermediate axis of the system. As expected, the 
Newtonian iso-potentials fall apart from the matter-iso- 
densities and acquire a form that is much closer to 
sphericity. In the scalar field cases, we find not only that 
their iso-surfaces follow the density distribution much 
closer than the Newtonian potential, but also that there 
are differences between them (even when the range of the 
field is the same for both models). 

We find that symmetron fields tend to be more spher- 
ical than the chameleons. The reason for such difference 
lies in the mechanism driving the screening of the fifth 
force: in the symmetron there is a defined threshold den- 
sity above which the field decouples from matter (the fifth 
force vanishes), while in the chameleon the fifth force 
disappears slowly and continuously (its range becomes 
shorter) as the density becomes higher and higher. 

To quantify the differences between the models and 
understand the dependence of the result with the model 
parameters, we calculate the axial ratios for every model. 
The values obtained for the Newtonian potential are 
q 2 = 0.84 and s 2 — 0.76, which differ considerably from 
the input parameters given for the density. The results 
obtained for the symmetron and chameleon models are 
shown in Fig 12] In order to make a fear comparison we 
show here the results as a function of the range of the 



field A s and A c that corresponds to the mean density of 
the universe in both cases. We find a complete differ- 
ent behavior for each scalar field model. The symmetron 
mechanism tends to give more spherical iso-surfaces when 
increasing both A s and zssb- The chameleon mechanism 
is insensitive to changes in A c up to ranges that are larger 
than the virial radius itself. Once again, this reflects the 
difference in the way the fields are screened. 

To test the stability of our calculations, we made res- 
olution and box size convergence tests by increasing the 
resolution by a factor of two while keeping the box size 
constant and also by increasing both box size and resolu- 
tion by the same factor. The tests were made for models 
parameters that are representative of iso-surfaces with 
small and large values of q 2 and s 2 . We find that our 
solutions are independent of the resolution. The New- 
tonian and symmetron solutions are stable with respect 
to changes in the box size. In the case of the chameleon 
model, we find that the solution is much more sensitive to 
the distance between the halo and the boundary. Never- 
theless, the actual change in the axial ratios when going 
from 16 to 32 Mpc/h in the box size is only of the or- 
der of 5%, which is far from the variations we see when 
changing from model to model. In any case, the environ- 
ment around non-isolated clusters is expected to change 
the solutions, but not the bulk of the signal. 

An estimation of the importance of the effect from 
the observational point of view can be made under the 
assumption that the X-ray component follows the iso- 
surfaces of the total potential (Newtonian plus scalar 
field). The ellipticity e = 1 — 6/aof the projection on 
the sky of these iso-surfaces should then be the same as 
of the gas density distribution. This assumption makes 
also possible to obtain constraints on the model param- 
eters, since there is a set of parameters for which the 
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FIG. 3: Contours of relative difference Ae/ejv between ellip- 
ticities that correspond to the modified models (LEFT: sym- 
metron, RIGHT: chameleon) with respect to Newtonian grav- 
ity. At small A the screening mechanism emerge and the al- 
lowed parameter space increases. Here zssb = 2.1 and n = 4. 



difference between the ellipticity that we predict by as- 
suming standard gravity and by including the fifth force 
is larger than the errors in present X-ray observations |9j . 

The calculation of the total potential was made by tak- 
ing into account that the geodesies equation at redshift 
z = has the following form in the symmetron case: 
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and can be written as: 
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for the chameleon. We estimate the ellipticity of these 
iso-surfaces of total potential by taking the mean value 
over random projections. Figj3] shows contours of the 
relative difference Ae/ejv between ellipticities that corre- 
spond to the total potential associated to the modified 
models and to the Newtonian potential. The regions from 
black to light grey correspond to values larger than three, 
two and one times cr b s /e b s . In other words, the black re- 
gion is ruled out with more than 3a b s /£ bs- When mak- 
ing this comparison, we assume that observations and 
Newtonian theory give the same values for the ellipticity. 
By taking into account a possible bias of the predictions 
of the standard model towards more spherical halos 
one can relax slightly the constraints and include models 
with higher values of (3 that are excluded here. 

From FigEl it is clear that for high couplings and small 
ranges, the chameleon model is less constrained than the 
symmetron. In other words, the chameleon tends to give 
more spherical objects. That seems to be in tension with 
Fig|2] in there, the symmetron field tends to be more 
spherical. This can be understood taking into account 
the dependence on zssb in the symmetron geodesies 
equation. This extra factor makes the total potential 
to have a stronger contribution from the scalar field and 



thus, a larger ellipticity than the chameleon, even in the 
case that its intrinsic distribution is more spherical. 

Notice, at galactic scales the iso-surfaces of the scalar 
fields should also follow the density distribution of the 
disks. Therefore, a signature of the screening mechanism 
may be measured in the MW via G AIA observations [l(| ■ 

In summary, we propose an astrophysical test which 
can be used as a probe to detect or differentiate screening 
mechanisms associated to scalar fields which are present 
in gravity theories which modified GR at scales larger 
than the solar system. We show that the existence of 
such screening mechanism can strongly affect the shape 
of galaxy clusters. Starting from a dark matter density 
distribution that corresponds to a non-spherical cluster 
of galaxies, we measure the shape of the iso-surfaces that 
correspond to Newtonian potential and two scalar field 
models (symmetron and chameleon). We find that both 
scalar field models give iso-surfaces that follow much 
more closely the density distribution than the Newto- 
nian potential. Furthermore, we find that the shape 
of the iso-surfaces also depends on the mechanism used 
to screen the fifth force: the symmetron model tends 
to give more spherical distributions than the chameleon 
one. Since, present observations show some tension be- 
tween the shape of real clusters 0, [ll[ and predictions 
obtained from simulations 12|, our results indicate that 
if scalar fields make any difference, it is in the right di- 
rection to correct the discrepancy in the observations. 
Finally, using recent data from X-ray observations, we 
calculate novel constraints on the coupling /3, and the 
ranges, A Si o and A c , of the scalar fields' fifth force. 

DFM and CLL thank funding from Research Council 
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